The consistency of codimension-2 braneworlds and their cosmology 
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We study axially symmetric codimension-2 cosmology for a distributional braneworld fueled by a 
localised four-dimensional perfect fluid, in a six-dimensional Lovelock theory. We argue that only the 
matching conditions (dubbed topological) where the extrinsic curvature on the brane has no jump 
describe a pure codimension-2 brane. If there is discontinuity in the extrinsic curvature on the brane, 
this induces inevitably codimension-1 distributional terms. We study these topological matching 
conditions, together with constraints from the bulk equations evaluated at the brane position, for 
two cases of regularisation of the codimension-2 defect. First, for an arbitrary smooth regularisation 
of the defect and second for a ring regularisation which has a cusp in the angular part of the metric. 
For a cosmological ansatz, we see that in the first case the coupled system is not closed and requires 
input from the bulk equations away from the brane. The relevant bulk function, which is a time- 
dependent angular deficit, describes the energy exchange between the brane and the six-dimensional 
bulk spacetime. On the other hand, for the ring regularisation case, the system is closed and there 
is no leakage of energy in the bulk. We demonstrate that the full set of matching conditions and 
field equations evaluated at the brane position are consistent, correcting some previous claim in the 
literature which used rather restrictive assumptions for the form of geometrical quantities close to 
the codimension-2 brane. We analyse the modified Friedmann equation and we see that there are 
certain corrections coming from the non-zero extrinsic curvature on the brane. We establish the 
presence of geometric self-acceleration and a possible curvature domination wedged in between the 
period of matter and self-acceleration eras as signatures of codimension-2 cosmology. 



I. INTRODUCTION 



Codimension-2 distributional defects 1 are infinitcsimally thin matter filaments which arc particularly interesting, 
subtle and elusive objects in gravitational theories. In four-dimensional general relativity, codimcnsion-2 defects of 
infinitesimal size describe cosmic strings which are the singular rclativistic versions of finite size vortices in condensed 
matter systems. Strings were thought to play a dominant role in the structure formation history of our Universe [1]. 
They were widely studied in differing contexts, but mostly at the zero thickness limit and in the test approximation; 
in other words, neglecting their self-gravity, i.e. the effect of their energy- momentum tensor on their proper evolution. 
When one takes into account their self-gravity, at least for a straight distributional string, one finds that they are 
described by a codimension-2 conical singularity [2] sourcing the constant tension of the defect (see also [3] for constant 
curvature bulk). Their effect is a global, topological one, since locally the induced gravitational field is flat. Indeed, 
each point of the string worldsheet generates a two dimensional cone with a non-trivial deficit angle. 

When considering however, non-trivial geometries, one was confronted by a paradox [4] which can be summed up in 
the following way: test cosmic strings are described using the Nambu effective action which leads to a local minimal 
area requirement, or equivalently a minimal surface spanned by the string worldsheet. This is geometrically described 
by the trace of the second fundamental form being equal to zero (for examples see [5]) and can be visualised as the 
surface (with boundary) of a soap film obtained upon retrieving slowly a circular wire from liquid soap. To the Nambu 
action one can find specific extrinsic and intrinsic curvature corrections [6] in order to take into account the finite 
width of the defect. It is natural therefore to expect, that upon considering the se/!f gravitating field of the cosmic 
string, a small correction should result on top of the minimal surface motion. Surprisingly however, self-gravity of the 
distributional string leads to trivial motion or non-bending of the worldsheet with the full second fundamental form 
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being identically zero. Geometrically the string worldshcet is a totally geodesic surface. In other words self-gravity 
makes the defect completely rigid not allowing it to bend in the ambient spacetime [7] giving back again essentially 
the straight string metric (see also [8]). 

This puzzling property of gravitational rigidity reemerged in the context of codimension-2 braneworlds embedded 
in a six-dimensional spacetime. Unlike four-dimensional General Relativity (GR), where one could convincingly argue 
that self-gravity of cosmic strings is negligible, this was far more subtle for a braneworld given that its motion is by 
construction fueled by its proper matter density. Indeed, not surprisingly, early work on codimension-2 distributional 
braneworlds showed that they were pure tension objects [9] . In other words, one could not find a distributional solution 
(see [10] for a mathematical discussion on distributional sources in GR) for a non constant energy-momentum tensor 
such as that, say, of a perfect fluid necessary to describe braneworld cosmology for example. One had to introduce 
finite thickness effects [11], which made the problem considerably more difficult [12], and more importantly, plagued 
the generality of the result and hence its physical relevance. Indeed, the key point when studying distributional 
sources is that they describe the important main features of brane dynamics for an arbitrary family of finite width 
regularisations at the limit of infinitesimal thickness. At the absence of a distributional description the dynamics are 
reguralisation dependent. The distributional description for example, leads to considerable simplification when one 
looks at codimension-one braneworlds where one can solve the full system of brane-bulk equations at least for cases 
with enhanced symmetry such as homogeneous and isotropic cosmology [13]. 

Activity on this subject concentrated mainly on the interesting topological properties of codimcnsion-2 defects in 
relation to the self-tuning paradigm [14]. The key point to resolve the gravitational rigidity puzzle was to understand 
that it was not the defect construction which was problematic, rather the gravity theory itself did not have the 
relevant differential complexity in order to describe complicated distributional solutions 2 . In other words, it is not 
strings that cannot bend, it is just Einstein equations that cannot describe the curving of strings infinitesimally. Even 
ordinary codimension-1 junction conditions fail in cases where the spacetime is of lesser symmetry, such as when 
pasting together the Kerr metric with flat spacetime, where the spacetime metric fails to be continuous. 

In fact, given that the defect carries a distributional energy-momentum tensor, the gravitational field equations have 
to be second order in such a way as to have piecewise continuous first metric derivatives. Although four-dimensional 
GR is the unique second order tensor theory with this property, in five or six dimensions one has to add an extra 
term, the Gauss-Bonnet term, to the action in order to have the most general second derivative (and not higher than 
second derivative) field equations. The generic second order derivative tensor gravity theory is in fact well-known to 
be Lovelock's theory [15] (for a review see [16], [17]). In this gravity theory each 2n + 1 dimensional manifold picks up 
in the gravitational action a novel Langrangian density, which is a specific combination of the nth power of spacetime 
curvature. This gravitational term stems from the Euler-Poincare characteristic of the manifold in 2n dimensions 
and hence is purely topological at this dimension. As such, the Einstein-Hilbert and the Gauss-Bonnet terms are the 
Euler-Poincare characteristics for two-dimensional and four-dimensional manifolds without boundary. 

The clear-cut hint that the completion of GR in higher dimensions may be relevant to codimension-2 braneworlds 
came with the work of Bostock et al. [18], where it was noticed that upon considering the general second derivative 
gravity theory, one could have, at least in principle, a non trivial energy momentum tensor fueling geometric junction 
conditions for a codimension-2 conical defect. Furthermore, it was shown that the higher order gravity term in question, 
the Gauss-Bonnet term, was generating on the brane an induced Einstein-Hilbert term plus extrinsic curvature 
corrections, whereas the Einstein-Hilbert term in the bulk was giving a pure induced cosmological constant. This 
observation was urged further in [19] (see also [20]) where it was shown that Lovelock theory can source defects 
with codimension even higher than two, unlike ordinary GR. A simple geometric explanation was given relating by a 
simple sum rule each nth bulk Lovelock density with the distributional brane's even codimension 2m and the induced 
brane Lovelock density, n — to. Therefore, the six-dimensional bulk Gauss-Bonnet term (n = 2) is "reduced" to a 
codimension-2 Einstein-Hilbert term on the brane n — m ~ 1 and so forth. The same n = 2 Gauss-Bonnet term in 
eight dimensions would " reduce" to a pure tension term on a codimension-4 brane n — m = and so on. 

This activity came to a halt when under some quite mild symmetry hypothesis (axial symmetry) it was claimed 
that the full set of junction plus bulk field equations at the location of the brane led to an inconsistent system of 
differential equations for a non-trivial distributional codimension-2 defect [21]. Furthermore, the junction conditions 
were shown not to be unique [19], without any obvious particular physical difference but with only differing mathe- 
matical regularity. Lastly, certain restrictive initial conditions on the braneworld surface [18] introduced important 
constraints [22] on the admissible matter on the brane. 

The aim of this paper is to falsify or explain all points in the previous paragraph. We will show that the system 
of equations is not only consistent, but has in general a free degree of freedom in the face of a varying deficit 
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angle function. This will result from carefully considering the correct and all relevant geometric expansions for the 
brancworld geometry. We will show that differing mathematical regularity leading to differing junction conditions 
inevitably leads to distinct physical setups (following up work by [20]). We shall show that the unique junction 
conditions tailored for a pure codimension-2 defect have the same mathematical regularity as those of GR and are 
the topological matching conditions of [19]. The matching conditions introduced in [18] inevitably lead to extra 
codimcnsion-1 distributional matter. Having established these basic facts, we will go on to show that one can obtain 
the modified Lemaitre-Friedmann-Robertson- Walker (LFRW) brane equations by solving the full system of the field 
equations, as well as the junction conditions in the brane neighborhood. We will find that the cosmological equations, 
apart from the ordinary LFRW term, have several interesting features including a geometric self-acceleration (in 
agreement with the maximally symmetric solutions of [23] ) , brane bending effects depending on the equation of state 
of matter, and energy exchange between the bulk and the brane triggered by a dynamical conical deficit angle. 

The organisation of the paper is as follows: we will give the general set-up and a careful treatment of the differing 
matching conditions (topological, topological with ring regularisation and general) in the next section. We will then 
demonstrate consistency of the whole set-up and go on to investigate the cosmological equations. We will conclude 
in the last section. 



II. GENERAL SETUP AND MATCHING CONDITIONS 

Let us consider the general system of six-dimensional Lovelock gravity coupled to localised brane sources. If we 
describe the regular bulk matter by S bu lk [g] and the distributional brane matter by S br ane [g] , the total action of the 
system is, 

S = d 6 x^\^\{n+^C GB )+S bulk [g] + S brane [g] , (2.1) 

with the Gauss-Bonnet Lagrangian density 

£gb =TZ 2 - AIZmn 1Z mn + TZmnka TZ mnka . (2.2) 

In the above, the calligraphic quantities refer to the bulk metric tensor g, while the regular ones to the brane metric 
tensor g. The field equations arising from the action (2.1) are 

pN a n-/N _ JlrrN , 2 rpN in n\ 

»M ~q^M 6 M ' 6 M ' y^- 3 ! 

where Tmn is a regular bulk energy-momentum tensor and Tmn is the distributional brane energy-momentum tensor. 
The most plausible assumption is that the bulk contains only cosmological constant Tmn — g MN , however, 

for the sake of generality we leave for the moment Tab arbitrary and non-vanishing (we will however assume later 
that it is regular at the position of the brane). The Gauss-Bonnet contribution to the equations of motion is explicitly 

Ti-MN = 2^gb Zmn^^TZ1Zmn+^TZmkTZ n k + AlZMKNATZ KA — 21ZMKA i zTZ N Kk " ■ (2-4) 

In the following, we will fix the notation to 3k§ = Air to simplify the equations. 

Let us now consider that there is axial symmetry in the bulk, so that the bulk metric ansatz can be written in the 
brane Gaussian-Normal coordinates as 

dsj = dr 2 + L 2 (x, r)d6 2 + g^(x, r)dx^dx v = dr 2 + g ab (x, r)dx a dx b , (2.5) 

where 9 has the standard periodicity 2ir, and the braneworld metric g^(x,0) is assumed to be regular everywhere, 
with the possible exception of isolated singular points on the brane. We denote by y collectively the transverse 
coordinates r, 9. 

To study in detail the matching conditions we must first differentiate between two distinct cases: 

1. The first case is topological matching conditions, discussed in [19], which have a geometric origin based on the 
distributional version of the Chern-Gauss-Bonnet theorem [24]. They assume everywhere smooth intrinsic and 
extrinsic tangential sections. 

2. The second case, which admits the lesser mathematical regularity, was introduced in [18] and assumes, not only 
a conical deficit based on the normal geometry, but also a distributional jump in the extrinsic tangential sector, 
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namely a specific combination of extrinsic curvature quantities 3 . We will call this possibility the general case 
since for regular extrinsic parts it reduces to the topological case outlined above. 

As we shall see later on, the second type of matching conditions inevitably lead to additional codimension-1 
distributional singularities [20], and hence additional matter sectors have to be introduced in order to close the 
system mathematically. This is not too surprising, since mathematical discontinuities in the extrinsic sector are 
expected to produce codimension-1 defects [4], [25]. On the other hand, the topological type of matching conditions 
lead to regular induced metrics unlike the former case. 



A. Topological matching conditions: Smooth regularisation 

Let us start with discussing the case of the topological matching conditions. In this case, there is no jump of the 
extrinsic curvature inside and outside the thin defect. As usually done when dealing with distributional sources, the 
matching conditions are derived by integrating around the singular space. Here, we seek to integrate over the conical 
space, say C, with internal metric 

ds\ = dr 2 + L 2 d6 2 , (2.6) 
with L(x, r) = 0(x)r + \h{x)r 2 + ^0 3 {x)r 3 + G(r 4 ) , (2.7) 

where 6 — 2n(l — f3(x)) is the conical deficit. Notice, the coefficients are allowed to depend on the "branc" coordinates 
x^. The angle 9 varies in the interval [0, 2ir). 

In the present paper, we are interested in the case that the distributional brane energy-momentum tensor supports 
a codimension-2 Dirac singularity and therefore can be written as 

T M N = T ^7^$M S N ■ ( 2 - 8 ) 

In general, there can be also codimcnsion-1 S(r) parts in the energy momentum tensor, which, as we will see later, is 
unavoidable if the extrinsic curvature has a jump across the interface layer of the defect. 

Then, by integrating the equations of motion, the only non-zero contribution yielding the codimcnsion-2 matching 
conditions reads [19], 



^ Jd 2 y^g- 2 R 2 



3 

la 



= T„ v , (2.9) 



where R 2 = —2L"/L is the curvature of the two-dimensional internal space, and a prime denotes d r . The Einstein 
equation is corrected with the following contribution from the extrinsic curvature K^ v = ^drg^ on the brane 

IV = K$K vX - KK^ + \g^{K 2 - K 2 kX ) , (2.10) 

which is identical to the one in [18]. 

In order to evaluate the integral of (2.9) which is singular at the tip, we assume a family of smooth regularisations 
C\ of the cone, parametrised by A, with smooth C 2 everywhere caps (for details sec [8]), pasted smoothly to the rest 
of the cone at a boundary dC, which we assume does not depend on the parameter A (see left part of Fig.l). The 
distributional limit is obtained at the A — > A c limit, where the cap tends to the singular conical space. The width 
of the regularised defect is supposed to be small enough so that the r-dependent brackets can be approximated with 
the leading term in their r-expansion, which is just the value of the brackets at r = 0. Then the integral has to be 
performed over only the internal curvature. The Chern-Gauss-Bonnct theorem relates the geometric curvature for 
arbitrary A, to the Euler characteristic \ 01 the regularised cap and the integral of the geodesic curvature k g of the 
boundary dC. Then, for every cap labeled by A, the following holds 



^- [ d 2 y^R x = X -^- I k g d6. (2.11) 



In [18] a particular jump of the extrinsic curvature was assumed, namely that the extrinsic curvature vanishes at the core of the defect. 
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FIG. 1: On the left, the conical singularity is regularised by a family of regular interiors C\ smoothly connected to the exterior 
solution of the cone at the boundary dC. On the right, a ring regularisation is chosen where the internal space function L is 
not smooth, but has a cusp across the ring at r = e. 



Since the A- labeled caps are topologically equivalent to a disk which has x — 1 and since k g = L', the theorem gives 
that for every kind of smooth regularisation of the conical singularity, one obtains 

/ d 2 y^R x = l-(3 . (2.12) 
4tt J Cx 

The above result only used the fact that the expansion of L close to the conical defect has the expansion (2.7) and is 
independent of the smooth regularisation of the interior of the defect. Finally, with the use of the above theorem, the 
codimension-2 matching conditions for the topological case is 

Let us now look at the leading 0(l/r) terms of the equations of motion. These come from the contributions of terms 
of the equations of motion which are multiplied by L' /L, or V^L'/i. With V M we denote covariant differentiation 
with respect to the metric g^ v . Then, according to the expansion (2.7), these terms contribute as 



L 



^f-l+Od). (2.14) 



With the natural assumption that the energy momentum tensor Tmn does not blow up close to the distributional 
singularity (otherwise the singularity would not be distributional), we obtain that the 0(1/ r) terms of the (rr) and 
the (r/i) equations yield respectively 

[g^ + VtV - ^<v) = > ( 2 - 15 ) 
" ' : ^ G^ u + Wjxu - ^-g,J) + = . (2.16) 



(3 V 2cT 

These two equations act as constraint equations for the two quantities K^ v , /?, which have to be determined and 
substituted back to the brane Einstein equation (2.13). In the general case, there are 11 independent unknown 
functions, but there are 5 independent constraint equations in (2.15), (2.16). Therefore, there will in general be 
functions which cannot be determined by the local equations around the brane, but need to be determined by the 
bulk solution. This, in fact, is clearly visible in the case of isotropic cosmology on the brane, as we will see in the 
next section. 

Finally, if we differentiate the brane Einstein equation (2.13) and use the constraint (2.16), we can derive the energy 
conservation equation 



From the above, we see that the energy is not strictly conserved on the brane, but can radiate in the bulk if the deficit 
angle changes. We can therefore anticipate that cosmological evolutions, close to the standard four dimensional one, 
will have /3 almost constant. We will come to this point in a later section. 
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B. Topological matching conditions: Ring regularisation 



In the previous analysis, we chose an arbitrary smooth regularisation of the distributional part of the geometry 
and the source. In this section, we will analyse a more specific case, where the codimension-2 source is taken by the 
limit of a sharp codimension-1 source with infinitesimal radius. Therefore, we replace the tip of the cone with a cap 
glued to the rest of the compactification with a ring interface (see right part of Fig.l). The extrinsic curvature is still 
continuous across the ring, but the angular part of the internal space metric L has a cusp. Inside the cap, since we 
will finally take the limit of its thickness to zero, the normal derivatives in the center of it tend to the values of the 
corresponding quantities just inside the ring. 

Outside the ring, we assume that the geometry has an expansion as 



oc ^ 



(2.18) 



n=2 



which is the standard conical singularity structure in the limit of zero ring thickness. On the other hand, inside the 
ring, the geometry should tend to two-dimensional flat space 



oo 

L-(x,r)=C(x,e)+r + J2^Mx)r n . 



(2.19) 



x 1 

with the constant C(x, e) = e[(3(x) — 1] + ^ — (f3 n (x) — /3 n (x))e n to ensure continuity of L across the ring. In the 

above, e is the width of the cap. As we will see later on, it is important for the consistency of the model that the 
higher order coefficients in the expansion of L (at least fa) arc different inside and outside the ring. The extrinsic 
curvature across the ring is continuous as in the previous subsection. We will assume that the derivative of the 
extrinsic curvature is also continuous, although this will turn up to be a fact, rather than an assumption, when we 
study the consistency of the setup. The crucial difference with the previous subsection is that the metric expansion 
inside the ring (2.19) is specified and L' has a definite jump. 

In general, there arc two types of distributional singularities that one can obtain in the e — > thin limit. One is 



S(r-e) 



There can exist, 



that of codimension-2 which has the structure <5^ 2 '(r) = and originates from the terms 
however, in the same limit a codimcnsion-1 defect with singularity structure S(r), which comes simply from S(r — e). 
To calculate the distributional pieces of the equations of motion, we will make use of the following identities where 
we denote by "distr" the distributional parts of the corresponding quantities 



e}6(r 



distr(i'V) = {P - l)[g^ + g'^ 
distr(L"G^) = (/3-l)[G^ 
distr(L"T4V) = (p-l)[W^ + W'e}S(r 



e) 



(2.20) 
(2.21) 
(2.22) 



Here, we also note the next to leading order terms in the e-expansion, since these distributional parts are divided 
by L ~ e at the ring position. Because of this well defined ring structure in the internal space, there could be 
codimension-1 contributions of the equations of motion surviving in the e — ■> limit. It is clear to see that the only 
such possible contributions come from the (/j,v) equations. To find the codimension-1 contributions, we should be 
careful with the distributional definition of Tmn in (2.8). This is because, inside the cap, the energy momentum 
tensor which becomes distributional in the e — > limit, may have an expansion of the type 



T^ = T^+T^r + 0(r 2 ) . 



Expanding (2.8) at r = e, we obtain 



-MN 



^^+(T«-i|T(0))^-e) 



1 

2^/3 



(2.23) 



(2.24) 



In the following we will use T^ v to denote T$ . 

From the 5(r — e)/e singular pieces of equations (2.3), we obtain the matching conditions for the codimension- 
2 singularity. Due to the topological origin of the dimensionally extended Lovelock densities, the Gauss-Bonnet 
combination in the bulk theory will give an Einstein brane term plus extrinsic curvature corrections, at the level of 
the junction conditions. The latter then read as before 



G 



i 



a(l 



T 

0) 



(2.25) 
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From the S(r — e) part of the (fxv) equations, we obtain the codimension-1 matching conditions at r = e 



2a 



a(l 



T (i) 

13) 



(2.26) 



As we will comment at the end, we need always have non-zero T$ in order not to overconstrain the system. Let us 
note that (2.26) is a consequence of the ring regularisation. Such kind of codimension-1 equation did not arise in the 
the smooth regularisation in the topological case, since then, only two-dimensional distributions are well defined. 

Let us now look at the 0(1/ e) terms of the equations of motion. We will make again the natural assumption 
that the non-distributional energy-momentum tensor Tab inside and outside the rings is regular and has no 0(1/ e) 
singularity. In evaluating the 0(1/ e) terms, we should be careful to do it both infinitesimally outside of the ring and 
also infinitesimally inside of it. This is because the function L' has a jump. This function enters in the equations 
of motion in two forms, as L'/L and as V^L'/L. Then, according to the expansions (2.18) and (2.19), these terms 
contribute as 



T 



T 



i 



i i 



L 



+ 0(1) 



e 

>—+0(e) . 



The (rr) equation has only L'/L pieces, therefore the 0(1/ e) terms give us the equation 

(g^ + - = . 



(2.27) 
(2.28) 

(2.29) 



On the other hand, the (r/i) equation has both L'/L and V^i'/i pieces. The difference of the contributions of these 
two pieces is a V"W M „ term, which therefore has to vanish 



vw^ - . 

The 0(1/ e) part of the the (rfi) equation just inside the ring gives by itself 



V/3 



(2.30) 



(2.31) 



Counting again the constraints of (2.29), (2.30), (2.31), we find that there are 9 independent equations to determine 
the 11 unknown quantities, and [3. Therefore, there will again be in general functions which cannot be determined 
by the local equations around the brane, but need to be determined by the bulk solution. However, in the case of 
isotropic cosmology on the brane, as we will see in the next section, all the unknown functions entering in the brane 
Einstein equation (2.25) are completely determined (modulo integration constants) and the system therefore closes. 

Had we taken T$ = in the codimension-1 matching conditions (2.26), which have 10 independent components, 
these matching conditions would act as constraints and therefore in principle overconstrain the system. This is indeed 
the case for the cosmological ansatz case that we will discuss later. 

Finally, let us again differentiate the brane Einstein equation (2.25). Using the constraint (2.30), we see that the 
brane energy momentum tensor is always conserved 



V7 F = . 



(2.32) 



Therefore, the ring regularisation, i.e. the fact that we consider not a C 2 but a C 1 cap C\ at r — e, freezes brane 
radiation. 



C. General matching conditions 

To study the general matching conditions, we will use a similar regularisation as in the previous subsection, where 
in addition, the tangential sector will have in principle distributional jumps. The regularisation of the codimension-2 
singularity is then depicted in Fig. 2. In more details, we assume that both the angular metric function L and the 
extrinsic curvature have jumps when going from the exterior to the interior of the ring. Inside the ring, since we 



FIG. 2: The conical singularity regularised by a cap glued to the bulk at a ring intersurface at r = e. Both and L have 
jumps across at the ring. 



will finally take the limit of its thickness to zero, the normal derivatives in the center of the cap tend to the values of 
the corresponding quantities just inside the ring. 

In this general setup, we would like to investigate the possibility of having as before a codimcnsion-2 distributional 
energy-momentum tensor as in (2.8). We consider the case of having a codimension-1 distributional source as out of 
the scope of this paper. The answer to the former question can easily be seen in the following to be negative. 

The extrinsic curvature inside and outside the ring brane have in general a non-trivial relation between them, 
depending on the details of the internal cap structure. Without loss of generality of our following analysis, we can 
make the simplified approximation that the two extrinsic curvatures are proportional as 



(2.33) 



with rj = const. In the following, we will simplify the notation denoting with the exterior extrinsic curvature K+ v 
(and similarly for other quantities constructed by it). Furthermore, we assume that higher orders in the r-expansion 
of the metric are continuous across the ring brane. This means in particular that the non-distributional part of K' 
is continuous, in correspondence with /?2 in the expansion of L. 

When looking at the distributional part of the equations of motion, since also K^ v can now have a jump, we should 
modify the distributional parts in (2.20)-(2.22) as following 



distr(LV)' = p-l)<?^ + (/?-^)<?>]<Kr-e) , 

distr(L'G^)' - [(p-^G^ + ip-riG'^Sir-e) , 

distr(L'WV)' = [(/3 - r, 2 )W^ + (/3 - l)W'^e]d(r - e) 

distr^ = (1 - rj)K^6(r - e) . 



(2.34) 
(2.35) 
(2.36) 
(2.37) 



The difference in the e-terms in the brackets has to do with the assumption that K'^ is continuous across the ring. 
All the derivatives in the above expressions are computed at the exterior of the ring. From the above we can first 
write down the codimension-2 matching conditions from the S(r — e)/e terms of the (jiv) equations of motion 



v 2 -P 

1-/3 



WW 



2a 



9V = 



(2.38) 



For t) = 0we recover the matching conditions of [18]. The codimension-1 S(r — e) part of the (^v) equations of motion 
provide us with the codimension-1 matching conditions which, in agreement with [20], give after some simplifications 



1 — {3 3 



+/3 
+/3 



7] — 1 



jJ.U 

V K V„L 



3/3 ?? 



2a r] -13 
UL 



1 (K„ v - Kg„„) -13 V 



77-/3 



K*W KV 



1 



77-/3 



L 



V^L PL _ Yj^ K + 



1 



G K \K K g^ v + R^ v + -RK^ — R KU K^ — R^ K K^ 



1 r (1 

a( V -f3) 



(2.39) 



Note that the above equation for rj = 1 yields the codimension-1 matching condition (2.26) of the previous section. The 
codimension-1 5(r — e) part of the (09) equation, on the other hand, will not have an energy-momentum contribution 
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because we assume that the defect is of pure codimension-2 nature. Therefore, in agreement with [20], we obtain 



( I V )K V 3 )K^W^ - ^(1 - rj)K ttv G llv = . (2.40) 



From the above constraint equation there are two possibilities. The obvious solution is rj = 1, and which gives us 
the topological matching conditions. There could be of course some solution for satisfying the above equation 
for rj ^ 0. However, as we will see in the following section for the cosmological ansatz, it turns out that these 
codimension-1 parts overconstrain the system. 

Regarding the 0(1/ e) terms of the equations of motion in this general conditions case, we obtain the constraint 
(2.29) from the (rr) equation outside the ring brane 

K"" (g» v + WV - ^2^) = , (2.41) 

however, inside the ring, because the extrinsic curvature has a jump, we obtain by substraction of (2.41) 

V ( V 2 - 1)K^W^ = . (2.42) 
Moreover, from the 0(1/ e) terms of the (rp) equation, we obtain inside the ring 

(g^ + V 2 W^ - JLs M „) - , (2.43) 

and for the difference of the corresponing equation outside the brane from the one inside it 

VW„„ = (r? 2 - l^j-W^ . (2.44) 

For general rj, we see that, apart from the constraints (2.41), (2.43) and (2.44) that we also had in the topological 
ring regularised case, there are 2 more constraints from (2.40) and (2.42). Thus, since the total number of constrains 
is 11, we can in principle determine the 11 unknowns, K^ v and (3. However, it turns out that in cases of symmetry 
as the cosmological ansatz that we will study in the following, the system for general 77 becomes overconstrained. 



III. COSMOLOGICAL EQUATIONS AND CONSISTENCY 

In this section, we will study the cosmological equations for a codimension-2 brane. We will first study in which 
cases the system of the brane equations studied in the previous section closes. We will then move on by checking the 
consistency of the next order in the e-expansion of the equations of motion. For this purpose we adopt the following 
metric ansatz for the four-dimensional metric describing LFRW cosmology 



ds 2 = -n 2 (t, r)dt 2 + a 2 (t, r)ds 2 K + dr 2 + L 2 (t, r)d6 2 , 
with the spatial metric for k = ±1,0 curvatures 

dp 2 



ds 2 



1 — Kp' 



+ p (d(j> 2 + shr (j> du 2 ) . 



(3.1) 



(3.2) 



The expansion of n and a is 



n(t, r) 
a(t, r) = a 



1 + Nr + -N 2 r 2 + -N 3 (x)r 3 + G(r 4 ) 
2 



1 + At + \a 2 t 2 + ^ 3 (x)r 3 + 0(r A ) 



(3.3) 
(3.4) 



where with a — a(t) we denote the scale factor on the brane. We will assume that the energy-momentum tensor on 
the brane is the one of perfect fluid = (— p, P, P, P). The energy density and pressure may contain contributions 
from the vacuum energy of the brane. We will not factor out these contributions until the next section. 
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A. The brane equations of motion 



Let us now write the system of brane equations for the different cases discussed in the previous section. By the 
term brane equations, we mean the matching conditions supplemented with the constraints from the 0(1/ e) parts of 
the (rr) and (r/j) equations, and the codimension-1 constraints when appropriate. 

For the case of topological boundary conditions with a smooth regularisation, we obtain that the equations (2.13), 
(2.15), (2.16) give the following system 



3a(l - 13) 
p + iP 



6a(l 

/ = 



0) 
N 

~A 



2AA + 2H(l - f)A 2 = 



H 2 
a 

- -+ 

a 

3P 
P 

$ 



a 
1 

2^ 



K 1 a2 

- + 2-*- A 



fA 2 



(3.5) 
(3.6) 
(3.7) 
(3.8) 



/3 3a(l-/3) ' 

In order to derive the equations (3.7) and (3.8), we have used the first one (3.5). The above system of equations has 
evidently one free function, which we can take to be (3(t). This function is expected to be fixed by the asymptotic 
dynamics of the bulk and cannot be determined by the local equations of motion. 

For the case of topological boundary conditions with ring-type regularisation, we find that (3.8) is split in two parts. 
The equations (2.30) and (2.31) now give 

f)A 2 = , 
0p = 0, 



2AA + 2H(1 



(3.9) 
(3.10) 



From the last one, we will in the following choose the solution $ = 0, i.e. that the deficit angle is constant, since the 
other possible solution demands that the brane energy density vanishes (p = 0). In this case, the system of (3.5), 
(3.6), (3.7), (3.9) and (3.10) is closed, the cosmology on the brane is uniquely determined, and is not depending on 
the bulk dynamics but only on the local structure of the equations of motion. 

The general boundary conditions on the other hand, have corresponding equations with the ones of the previous 
case, plus two more which have to do with the constraint (2.42) and the codimension-1 constraint (2.40). These read 

P I-/? 



3a(r/ 2 - (3) 

P + 3P 
' Qa{n 2 - (3) 
N 
1 



1 



(H 2 + 



K 1 

? + 2^> 



+ )-fA 2 
rj z — p a 2a 

3P-3aA 2 (n 2 - 1) 

p+QaA 2 {rj 2 - 1) 



2AA + 2H(l- f)A 2 

/3( P + 3«/3(l -/3)(t7 2 - l)^ 2 ) 
r ? (r ? 2 -l)(l + 2/) 



(l-r?) 



2a 



f + 3 - — (1 + n + il 2 )A 2 (l + 2f) + 2aH 2 {\ + /) - 4a- 



t 






0. 



( V 2 -1)A 2 



A 2 (3.11) 
(3.12) 
(3.13) 

(3.14) 

(3.15) 
(3.16) 

(3.17) 



Obviously, for the topological case we arrive at the generic conclusions we explained in the previous section. Let us 
concentrate on the general case. From (3.16) we see that can either have r\ — or / = — 1/2 with general r\. The 
choice rj = is inconsistent for general matter on the brane because of the codimension-1 constraint (3.17). On the 
other hand the choice / = —1/2 fixes the equation of state to w = —1/6 (for tensionless branes - for branes with 
tension the equation of state will be fixed to another constant value), therefore it is again not so interesting. 

From the above we deduce that the only consistent boundary conditions with a purely codimension 2 defect are 
the topological ones. This theorem lies on two assumptions: first that the cosmological ansatz (3.1) does not contain 
any g tr or g r0 components (see [26] for criticism of the (3.1) ansatz) and second that the metric is analytic at r = 
so that it can be expanded in powers or r as in (2.18), (2.19), (3.3), (3.4). Let us also stress that, if we had had a 
true codimension-1 defect with non-zero Tgg, there would be no consistency issue for any ij, since then the constraint 
(3.17) would have a matter part in the right hand side. 
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FIG. 3: The grouping of equations according to the various functions that are to be determined. On the left table, the 
first column has the codimension-2 junction conditions, while the other columns have the various orders in r of the Einstein 
equations. The 0(1/ e) component of the (99) Einstein equation vanishes identically. 



B. Consistency of the next order equations 

Up to this point, we have looked at the matching conditions supplemented by some equations of order 1/e. These 
included apart from the brane metric, the deficit angle [3 and the extrinsic curvature K^ v . In the cosmological case 
the corresponding functions were the following: a, (3, A and N. It is easy to see that in general the various terms in 
the e-expansion of the equations of motion break into groups according to the unknown functions that they involve. 
The resulting grouping can be easily seen from Fig. 3. 

The next group of equations (which we call second order) is the one with the 0(l/e) terms of the (pv) equations 
and the 0(1) equations of the (90), (rr) and (rt). These form a set of five equations which involve the next order 
variables A 2 , N 2 and [3 2 . We will study them in the cosmological case to see if they are consistent. The strategy is 
to solve for them from three of these equations and then check that the two extra equations do not overconstrain the 
system. This check is done in view of [21], where it was claimed that the system is inconsistent (overconstraincd) . 
Here, we find that a careful inclusion of the (V ' /VjK^ terms omitted in [21] gives consistency. Without loss of 
generality, we will make the simplifying assumption for this proof that the bulk is empty, i.e. Tmn — 0, although 
some matter could be necessary in the interior of the defect for the ring regularised case. We will comment at the end 
of the section about what happens if we drop this assumption. 

Let us look first for the topological case without a specified regularisation. We can first solve for A 2 , N 2 and (3 2 from 
the (D(l/e) parts of (tt) and (ij) equations and from the 0(1) part of (99) equation. As it can be seen in Appendix 
A, these equations of motion form a 3 x 3 system of linear algebraic in the second order variables equations. Using 
the brane equations (3.5)-(3.8), we can easily solve for A 2l N 2 and (3 2 as 

15(1 - /?) 2 + 2(1 - (3)(p + 9P) + ApP - 

M 4a(l-/J)(„ + 9P) (3 ' 18) 

-p(p + 3P) 2 - f (1 - Pf(p + 6P) + 3(1 - f3)p(p + 9P) - 6p 2 P + 

N2 = 6a(l-(3)p(p + 9P) (3 - 19) 

45(l-/3) 2 -2p(p + 3P)-^£ 

A = ^ 2MP + 9P) ■ (3 - 20) 

The next step is to make sure that the 0(1) parts of (rr) and (rt) equations do not overconstrain the system. These 
equations are written down in Appendix A. It is a tedious, but straightforward exercise to see that in fact they are 
automatically satisfied and therefore the system is consistent. 

The next check we have to do concerns the ring regularisation in the topological matching conditions. Here we need 
to compute the equations also in the interior in order to see if the system is consistent. In this case, a simplification 
which will help us is that the equations for the first order variables imply that $ = 0. This simplifies considerably 
the equations. In the interior, having assumed that A 2 and N 2 are the same (continuity of K' ), we have the extra 
variable /3 2 . The algebraic system in the interior of the 0(l/e) parts of (tt) and (ij) equations and from the 0(1) 
part of (99) equation, is the same as for the exterior with the only substitution 

h. = j ■ (3.21) 



The next step for the consistency of the system in the interior is to see the 0(1) parts of (rr) and (rt) equations. 
The 0(1) part of the (rt) equation is indeed then identically satisfied. However, the 0(1) part of (rr) equation is not 
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satisfied without the inclusion of a 7^ r m ' ) matter in the interior. In more detail, the matter we have to include in the 
interior of the ring in order not to overconstrain the system is given by the latter equation to be 

r( in) _ -45(1 - (3f + 6(1 -(3)(p- 3P) - 2p(p + 3P) 

T - " = (3 - 22) 

Had we not taken a Trr^ , then the equations would only be satisfied for a specific relation of p and P of the brane 
matter. This would lead to overconstraining the system, and in the spirit of [21] to inconsistency. However, it is 
natural to assume that the thickening out of the thin codimension-2 defect gives rise to a non-zero energy momentum 
in the interior of it. 

As far as the third order equations of the table in Fig. 3 are concerned, there is no question about consistency. This 
is because there are three new variables A3, N3 and (3% with two equations (in the smooth regularisation case), the 
O(l) parts of the (tt) and the (ij) equations. In the ring regularisation case, there is one more variable /?3, but two 
more equations, as the number of equations is doubled (evaluated inside and outside the ring), so again every function 
can be determined. Therefore, we will not study them in more details. 

Finally, let us comment about the differences in the above approach if we had allowed for general matter T MN in 
the bulk. It is obvious that in this case the above solutions for the second order variables are altered. Then, for the 
smooth regularisation, we would obtain from the 0(1) parts of the (rr) and (rt) equations that 

T r8 = , 77 = % e . (3.23) 

The same would hold for the exterior to the ring matter in the ring regularisation case. For the interior, the equations 
of motion are satisfied if 

T r f l) = , T e {m) 9 = 1* , r = X + hi* , (3.24) 

where X is the interior rr component of the energy-momentum tensor in the absence of other matter given in (3.22). 



C. Vanishing extrinsic curvature limit 

In the previous section, one could wonder about what happens in the limit of vanishing extrinsic curvature (K^), 
since several expressions are divided by A or N. In this subsection, we will look back to the original equations 
without dividing them by components of the extrinsic curvature. First the brane equations read for both types of 
regularisation that we considered 

P tt2 , K , 1 



H z + — + — , (3.25) 



3a(l - (3) a 2 2a 

P + 3P (3.26) 



6a(l -[3) a 2a 

$ p = . (3.27) 

From the above we have two distinct choices. First that p = 0. Differentiating the Friedmann equation and using 
the acceleration equation, we find then that also P = 0. The solution then for the scale factor is de-Sitter space in 
all possible foliations (n = ±1,0). The second order constraints for the 0(1) parts of the (00), (ij) equations and 
the 0(1) part of the (rt) equation are trivially satisfied. On the other hand, the 0(1) parts of the (88) and the (rr) 
equations are only satisfied with the inclusion of bulk matter. In both types of regularisations (smooth and ring) and 
for both inside and outside the ring in the ring regularisation case, the later equations give 

T; = r; = il, (3.28) 

which corresponds (when also looking at the 0(1) parts of the (pv) equations) to cosmological constant k§ = —h/2a. 
This value of the bulk cosmological constant in fact realises the Born-Infeld limit [27] as noted in [23]. 

The second distinct case is when = 0. Notice that then, the equations (3.25), (3.26) are the 4-dimensional 
standard LFRW equations 4 . Then, the second order constraints for the 0(1) parts of the (00), (ij) equations and the 



4 Parts of these observations can be found in the PhD thesis of Paul Bostock [28]. 
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0(1) part of the (rt) equation are trivially satisfied for 2 — 02 = 0. On the other hand, the 0(1) parts of the (66) 
and the (rr) equations give constraints. In the smooth regularisation, they coincide and give a relation between A 2 
and N 2 . In the ring regularisation, the same holds for the outside problem. Inside the ring, the 0(1) part of the (66) 
equation gives the same relation between A 2 and N 2 - On the other hand, the 0(1) part of the (rr) equation needs 
some extra T rr . 

It is instructive to see the limit of this special case for Minkowski branc vacua, where H = and k = 0, with 
cosmological constant Aq in the bulk. Then, the brane has tension T4 = p = —P = 3(1 — 0)/2. In the smooth 
regularisation case, the relation between A 2 and N 2 is rather simple 

N 2 = -3A 2 - A 6 . (3.29) 

In the ring regularisation, the 0(1) part of the (rr) equation dictates that 

(3.30) 

where the first addendum is the bulk cosmological constant contribution and the second is an extra contribution 
necessary for the regularisation to work. We can compare then, this special case with the known exact solution 
of [23] of the double Wick rotated black hole. As we see in Appendix B, the above reproduce the correct relation 
between the bulk cosmological constant and the coefficients of K'. This example is important, because it violates the 
assumption of the form of the expansion of K^ v that was considered in [18] (namely K'^ = 0). The same assumption 
was subsequently used in [22] and led to erroneous constraints for the brane matter. 



q-(in) _ 



47T 



A fi + A, 



1-/3 



IV. CODIMENSION-2 COSMOLOGICAL EVOLUTION 



Let us now revisit the brane equations and try to understand what kind of cosmological evolution can be obtained 
on the codimension-2 brane. Let us stress here that, although we proved that the equations of motion locally around 
the purely codimension-2 defect are consistent for the topological boundary conditions, one has to make sure that 
there are no singularities introduced when integrating the equations of motion away from the defect. Here, we will 
not carry out this formidable task of the bulk integration, but we will content ourselves with deriving the possible 
four dimensional cosmologies from just the local equations on the brane. As discussed in the previous sections, the 
brane equations of motion for the smoothly regularised topological case are 

if 2 + 4 + 7 l-A 2 , (4.1) 



3a(l - 0) a 2 2a 

P + 3F 4 'i-/AV (4.2) 



6a(l - 0) a 2a 

2AA + 2H(l - f)A 2 = § — ^ -. (4.4) 

v 3a(l-0) K ' 

Here, we recognize the first two equations as a modified version of the four-dimensional Friedmann and acceleration 
equations. Note that is generically a function of time which can make the effective four dimensional Planck scale 
time varying. A second important point is the appearance of an effective cosmological constant in the face of a, which 
is of geometric origin. A last difference from the standard four dimensional equations is the presence of the extrinsic 
curvature correction parametrised by A. 

On the other hand, the continuity equation (2.17), which is not independent from the above, gives 

P + 3F(p + P) + ^=0. (4.5) 

Note that, in general, there is not energy conservation for the brane matter since a varying means inevitable energy 
exchange between the bulk and branc. This is in contrast with the 77 = case where there is energy conservation [29]. 
The energy conservation equation (4.5) can be rewritten as 

*( PM , 3 tf^±^=0, (4.6) 



dt\l-0J 1-/3 
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which looks like the standard energy conservation equation for a redefined energy density and pressure as (p, P) — > 
(p, P)(3/(l — (3). In fact, the five (four independent) coupled equations of motion, (4.1), (4.2), (4.3), (4.4), (4.5), which 
are the full information available to us at r = 0, are consistent with each other and constitute a non-closed system. In 
other words, one needs extra information coming from the bulk geometry in order to fix one of the functions, e.g. (3, 
and then to solve fully the system. The solution in the bulk is no longer unique as in the case of codimension-1 branc 
cosmology and one has a family of bulk solutions parametrised by the angular deficit function (3. Not all these bulk 
solutions will be acceptable. Certain of them will inevitably carry singularities away from the brane for example, and 
not only in the position of the brane at r = 0. 

From the matching conditions (2.13), we can read the effective four-dimensional gravitational constant as 

< " (£) M^W) • (47 » 

where we have momentarily reintroduced the six-dimensional gravitational constant according to (2.1). As noted 
before, an energy exchange between bulk and brane means a time- varying gravitational constant through the variation 
of (3. This variation is constrained during the early cosmology by the primordial abundances at the nucleosynthesis 

epoch, with a limit of this variation approximately < 0.2 for G = %-kk 2 (see [30] for details of this limit). This is 
constraining the variation of (3 as 







(1 - (3)H 



£0.2, (4.8) 



which is not a rather strong constraint. Further constraints come from the fact that the theory with varying (3 is 
similar to a scalar-tensor theory and therefore there will be strong constraints from solar system observations. Not 
knowing the full family of solutions in the bulk, we choose to consider the case where (3 is approximately constant. 
This, in fact, will be exactly the case if we restrict our analysis to the ring topological regularisation, where there is no 
energy exchange between the defect and the bulk. This case is bound to give us, at least seemingly, acceptable four- 
dimensional cosmology, and has the merit that the system of equations is closed and does not depend on undetermined 
(by the matching conditions) functions. 

When (3 is constant, we can solve for the extrinsic curvature A for a given, constant equation of state for the brane 
matter. For that purpose, we split the energy density and pressure to a part which has the form of brane vacuum 
energy A and part which describes brane matter 

P = A + p m , P = -A + P m , (4.9) 

with the brane matter equation of state defined as w — P m /pm- Then, for w ^ —1, we can integrate (4.4) as 

( " pW^ 1 , (4.10) 



2 



(A + p m ) 

with C 2 a positive integration constant. With this expression, the Friedmann and acceleration equations (4.1), (4.2) 
become 



K f 1 \ C 









) + 


2a 




K 2 


*( 


K 4 \ 




3 



p^ , (4.11) 



1 \ „ c 



^2 



3t^ ^(wp m ~ \)p% 1+w) , (4.12) 



a 6 V \i 2a J (A + p m f 

which show a non-trivial correction to the cosmological equations as a function of w. The magnitude of this correction 
depends on the integration constant C 2 . In the following, we will try to sec whether it is possible for large enough 
C 2 , consistent with observations, to have some interesting modification to cosmology. From the integration of the 
continuity equation (4.6) for (3 =const. we obtain 

which holds also for the special w = — 1 case. 

For C 2 7^ 0, our six-dimensional bulk geometry has a genuine curvature singularity at r — (apart from the 
distributional one). In fact, this is to be expected from purely geometrical considerations (see [31]). Higher codi- 
mension defects, when considered in their zero width limit, will develop curvature singularities. These are expected 
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FIG. 4: The functional dependence on the scale factor of the energy density and the extrinsic curvature correction to the 
Friedmann equation for different equations of state of the brane matter in the self-accelerating case. From top to bottom we 
have the equations of state for cosmological constant, domain walls, cosmic strings, dust, radiation and stiff matter. 



to be smoothable once we take finite width corrections into account. Therefore, we expect finite width effects to be 
important at the UV sector and our distributional approximation to break down even though this will not show up 
necessarily in the field equations themselves! 

Note that equations (4.11), (4.12) are valid not only for (3 exactly constant, but also for a (3(t) which is slightly 
oscillating around a constant value. To see this, we note from (4.6) that (as far as p rn + P rn > 0) the quantity 
u = p/3/(l — j3) has opposite monotonicity than a, and therefore, u can be used as a good time parameter. Note 
that for this case the solution (4.13) is not valid in general. From equations (4.1)-(4.4), we can obtain a differential 
equation for d(H 2 ) / du containing also the function j3. Using the small oscillatory behaviour of /3, we get an autonomous 
equation for d(H 2 )/du, whose integration gives equations (4.11), (4.12). 

Let us now analyse some particular cases of interest. 



A. Self-accelerating branes 

The first interesting case is the one in which the tension of the brane is vanishing A = 0. In this case, we see from 
(4.11) that there is an asymptotic accelerating phase for a < 0. This acceleration is purely due to the geometric Gauss- 
Bonnet term and therefore is a case of self-acceleration. This corresponds to generalisation of the self-accelerating 
solutions of [23] with the addition of matter. Before proceeding, let us note that in order that these self-accelerating 
solutions are viable, one has to check that they are free from ghosts that are typical in five dimensional self-accelerating 
braneworlds [32]. The correction to the Friedmann equation is given by 

A 2 = C 2 P I^ = J-^ . (4.14) 

Note that the second equality holds also for w = — 1. In the table in Fig. 4 we list the form of corrections that we 
obtain for different constant equations of state. It is worth noting four interesting limits of the cosmological evolution 
in this case. 

First is the one where the brane is dominated by cosmological constant. This may go against the initial assumption 
that the tension of the brane vanishes, because of the ambiguity to separate the vacuum part of the energy density 
from the matter energy density. If a matter component behaving as inflaton or dark energy dominates at some period 
of the history of the Universe, it will have the behaviour that we note here. So, for this case the extrinsic curvature 
plays the role of a dark radiation squared term a~ 8 that will dominate the cosmology at early times. 

A second interesting limit is when the extrinsic curvature on the brane vanishes ^4 = 0. This limit is of particular 
interest since we have a regular six-dimensional geometry at the location of the brane at r = 0. Then, the codimension- 
2 cosmology is exactly the four-dimensional LFRW plus an effective geometric cosmological constant which is 
positive for negative a. This is the generalised version of the solutions of [23] which were obtained in the case of pure 
tension. We see clearly that even in the absence of matter p = and n = we have a non-zero H 2 , hence these 
solutions are self-accelerating. 

A third interesting limit is when the matter equation of state is that of radiation w = 1/3. Then, the extrinsic 
curvature correction is of the form of constant vacuum energy. One would be tempted to use this vacuum energy to 
drive early Universe inflation. This seems to be possible, but around matter-radiation equality the Universe will be 
dominated by this new vacuum energy, making the phenomenology of the model problematic. 

A fourth interesting case is when the matter equation of state is that of dust w = 0. Then the extrinsic curvature 
correction is of the form of curvature (a~ 2 ). In this case, we can see a possibility of having an observable signature 
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FIG. 5: The log-log evolution of the energy densities as a function of the scale factor. With solid lines, there are the contributions 
of the standard matter and the geometrically induced 1/(2|q|) cosmological constant. With dashed lines are the evolutions of 
the extra component (oc A 2 ) in the Friedmann equation due to the extrinsic curvature of the brane. The dependence of this 
extra component on the scale factor is noted. This dependence changes whenever a new matter component becomes dominant, 
i.e. at matter domination at a eq and at the cosmological constant domination at a acc . 



of the codimension-2 cosmology, if the constant C 2 is chosen so that there is a brief period of curvature domination 
around the vacuum energy domination period. This scenario is depicted in Fig. 5. In order to obtain the standard 
epochs, i.e. matter domination and then cosmological constant domination, we need to have that p m {o, eq ) > l/(2|a|) 
and (K 2 /3)p m (a eq ) > A 2 (a eq ), where a eq is the scale factor at radiation-matter equality and a acc the one at the 
moment of cosmological constant domination. In order that the brief period of curvature domination is observable, 
we need A 2 (a aC c) > l/(2|a|). Therefore combining these inequalities we see that we need 

1 k 2 
^j- < A 2 (a acc ) < yp m (a e ,) . (4.15) 

Then, a straightforward manipulation with the help of (4.13), (4.14) provides us with the following constraint on C 2 

Po(W?A ~ 2/3 a 2 acc 2 k\ ( Pa (l~f3) \ 1/3 / a 2 acc \ 

2H < C < T l^lj • (4 - 16) 

It is not yet clear if such a brief period of curvature domination is phenomenologically viable and it requires certainly 
further study to see the potential observational signatures of such a case. 



B. Self-tuning branes 



A second interesting case is when self-tuning can be realised. For that possibility, the topological quantity (3 is 
used to accurately cancel the vacuum energy contribution of the brane to the effective cosmological constant. In more 
details, as discussed in [23], by tuning /? one can make the two vacuum energy contributions, the geometrical one 
1/2 1 a | and the brane tension one A, cancel 

This limit is orthogonal to the self-accelerating case that we studied before. Such kind of self-tuning solutions, found 
in [23] , have of course to be checked for their stability. The effective Friedmann equation for the present case is given 
by the expression 

„2 



c 2 



a? = 3 Pm + ( ^ +0 * < ( 4 - 18 ) 
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which includes a non-trivial correction beyond the standard linear to energy density term. Let us note here that the 
actual self-tuning is not visible from this Friedmann equation. This is because the present is valid for (3 = const. 
Instead, one should consider the time- varying deficit angle case in order to sec how (4.17) can be dynamically achieved 
at the early Universe evolution. Since before nucleosynthesis there are no constraints on the variation of Newton's 
constant, one could have an acceptable cosmology with varying (3. However, then, the brane equations do not close 
and the mechanism of self-tuning should be dictated by bulk boundary conditions. The study of such a case is beyond 
the scope of the present paper. 

V. CONCLUSIONS 

In this paper we have demonstrated that a distributional treatment of codimcnsion-2 branes is consistent and possi- 
ble if one considers the full Lovelock gravity in six dimensions. Higher dimensional GR cannot describe infinitesimally 
thin codimension-2 defects 5 for it does not present the differential complexity to allow such distributional terms [19]. 
Furthermore, we have shown that it is the higher order gravity terms in the bulk action that guarantee the presence 
of ordinary GR gravity on the brane. This, up to now, was hinted by the general form of the junction conditions [18], 
[19], but never demonstrated as a consistent solution of the full bulk and brane field equations. A nice example where 
the importance of the full bulk field equations becomes essential is the "apparent" choice of junction conditions [19]. 
In principle, at the level of distributions, differing mathematical regularity can give differing matching conditions, 
both a priori completely consistent. It is only upon considering the bulk field equations (in the spirit of [20]) that 
we see that the general junction conditions introduced in [18] inevitably need codimcnsion-1 matter sources. We 
have therefore established a genuine physical difference between the topological and general junction conditions. The 
former are pure codimension-2 junction conditions, whereas the latter are mixed codimension-1 and codimension-2 
junction conditions. It would be interesting to see the relation of the general matching conditions with work on 
intersecting branes and cascading cosmologies [33]. 

We have studied the cosmology of conical branes and have found several interesting features. Firstly, quite gener- 
ically the higher order terms give ordinary LFRW equations, as though one included an induced gravity term on 
the brane [34]. This is quite different from codimension-1 cosmology [13], where the ordinary behaviour is recov- 
ered only at late times. Corrections to this standard evolution are threefold: there is a geometric acceleration scale 
related to 1/a, which restates the problem of a minute cosmological constant as a gravitational hierarchy between 
the bulk Gauss-Bonnet and Einstein-Hilbert term. Secondly, extrinsic curvature corrections are apparent and they 
are dependent of the equation of state of the perfect-fluid matter. In other words, as the brane evolves in the bulk, 
a dust equation of state can lead to an extra component behaving as curvature, or a radiation equation of state to 
a cosmological constant term, and so on. To put it in a nutshell, each matter fluid introduces two differing fluid 
components in the modified LFRW equations. Last, but not least, as the brane evolves it can generically radiate in 
the bulk. Again, this is unlike the vacuum bulk codimension-1 braneworlds. 

The above characteristics can be used to constrain the model in question with respect to cosmological observations. 
Furthermore, on the theoretical side it would be interesting to have particular bulk solutions manifesting the cosmology 
evolution we have found, and in particular setting the boundary conditions in order to fix (3. The most important 
information we have given here in this direction is that it is pure six-dimensional Gauss-Bonnet gravity that gives 
the dynamical LFRW-like evolution on the brane. The six-dimensional Einstein term gives the possibility for a 
phase of geometric acceleration. Furthermore, a varying j3 could have interesting consequences for a cosmological 
self-tuning scenario following equation (4.17). However, note that the evolution in (3 would then leave an imprint on 
the cosmological evolution equations. These are amongst the open, interesting questions that this work puts forward 
and which we hope will be answered in the near future. 
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APPENDIX A: THE SECOND ORDER EQUATIONS OF MOTION 



In this Appendix we will list the equations of motion which involve the second order variables A 2 , N 2 and (3 2 . As 
we saw in the main text, they depend on the matching conditions and also the regulatisation that we have taken. For 
this Appendix, we assume that the bulk is empty, i.e. Tmn = 0, and furthermore that the three dimensional space is 
fiat, i.e. k = 0. 

Let us first see for the topological matching conditions without an explicit regularisation. The 0(1/ e) parts of (tt) 
and (ij) equations and the 0(1) part of (69) equation form an algebraic set of equations for the three second order 
variables. We can write the system in the suggestive form 




C 



(Al) 



with the following coefficient matrix 

( 2 

D= 2(1 + A) 

\ 2 -r{H 2 + ^-A 2 -2AN + 2-i) 
and the following vector of the right hand side 
/ 



C 



H 2 



A 2 + Z-AN+^ 

a 1 3 



' N 



2a 



(H 2 + h A 2 ) 



H 2 + ±-A 2 



H 2 ' ' 
A 3 N + H 2 (2A 2 



^-3A 2 + 24(^ 

2a N V a 



^s(H 2 + ^-A 2 



A/3 

W- p (H 2 + ±-A 2 -2AN + 2-i) 




AS 

2a I 



2N 2 



SAN) + 4(A + N)AH + 2A 2 + ^(A 2 + H 2 ) 



J 



(A2) 



(A3) 



The solution of this system is given in Sec. Ill B. There are two more equations that have to be verified from the 
grouping of Fig. 3. First, the 0(1) part of the (rr) equation. This is again an algebraic equation in the set of 
(P2,N 2 ,A 2 ) variables 



1 - 2aA(A + 2aN) 



2a 



H 2 



1 



H 2 



1(3 



A 2 +[l- 6aA 2 + 2aH 2 ] N 2 



+ 



A 

3P 



- - olA 2 - 3aAN H + aH 2 + a—H 2 + 2a- 

2 2 A A a 



P2 = 



AN + -N +- + ---- 2aA 2 A 2 + -AN + N ) + 2aH z -A 2 + -AN + N 



a 3/3 



--in a ( aA + ^nA- ^an 



i 



i 



2a d , , o 
3 a 



AAN + H 



la ■■ 2a /3 2 2 

y aa+ y/ h - A ) 



(A4) 



On the other hand, the 0(1) part of the (rt) equations is a differential (Bianchi) equation with respect to the second 
order variables 



4aA 2 - 



2a 
~Afl 



H 2 + — -A 2 
2a 



f3 2 +4a 



A 
A 



N 
~A 



H\A 2 - 4aHN 2 + 



2aA 



H I 



N 
~A 



A 



P2 - 



3H (l - + 3^ + 2aH [h 2 (l - + 3A 2 - 2N 2 - AN^j + 2a^(H 2 + 3A 2 + 2AN) . (A5) 
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For the case of ring regularisation of the topological case, we have as second order variables the set (A2, N2, 02, 
02). In the exterior of the ring, the system of equations in the same as above, with the simplification that = 0. In 
the interior of the ring, by taking the simple relation (3.21), we can trade 02 for 02- Then, we can easily see that we 
have again the same system (Al), with (A2), (A3). On the other hand, the 0(1) part of the (rr) equation changes to 



1 - 2aA(A + 2aN) 



2a 



H 2 



A 2 + [1 - 6aA 2 + 2aH 2 ] N 2 



A 



- - aA 2 - 3aAN + — - + aH 2 + a—H 2 + 2a- 
2 2 A A a 



0H - 0AN + —N + (1 - 0)A Z +0-- 2aA z ( A + 

Qj 



(6-/3) 



AN + N 2 ) + 2aH 2 [ I A 2 + 



(2-/3) 



AN + N' 



1 



3 3 

Finally, the 0(1) part of the (ri) equation also changes to 



AaH ( AA + ^NA - ^AN ) + ^- ((2 - 0)A 2 + AAN + 0H 2 ) - ^-AA . (A6) 



4a • 2a 
T A0 2 



H 2 + — -A 2 
2a 



02 + 



4a / A 

T 



+ 2- 



N 



H\A2- T HN2 + 1 p 



H[l + 



N 



3H 1 



N 
~A 



• 2,^ll \ II- ( J 



N 
1 

A 



! \A 2 - ^N 2 



1 AN 



+2a-(H 2 + [ — — l\A + -AN) . 











(A7) 



APPENDIX B: COMPARISON WITH THE DOUBLE WICK ROTATED BLACK HOLE SOLUTION 

In this Appendix we will compare the explicit simple solution of a system where the bulk has a Gauss-Bonnet term 
and cosmological constant, with the constraint that we noted in sec. IIIC. We will show that firstly, in this solution 
it is K' 7^ 0, violating the assumptions of [18] and furthermore that the leading terms in the r expansion of K^ v is 
the one found in the main text from the local consistency conditions. 

The solution of a Gauss-Bonnet AdS black hole with horizon of toroidal topology is 

ds 2 = -V(R)dt 2 + ^— + R 2 5 mn dx m dx n , (Bl) 
V{H) 



1 



\A- to (* 2 -£) 



(B2) 



with k 2 = — A 6 /10. Let us assume that a > 0, \i > 0, 4afc 2 < 1. Then, there is only one singularity at r = 0, shielded 
by an horizon at Rh = {^/k 2 ) 1 ^ . 

Let us now make a double Wick rotation t — > i0 and x° — > it. Then, we obtain the metric 

ds 2 = R 2 r ll _ lu dx"dx 1 ' + + V(R)d9 2 . (B3) 

V(R) 

We can go to a Gaussian-Normal system by means of the transformation R' = W and obtain the metric 

ds 2 = R{r) 2 r ltlu dx^dx u + dr 2 + V(R(r))d6 2 . (B4) 

Then, the comparison with the previous presentation is straightforward 

L = W , (B5) 
Kfjtu = \{R 2 )'ri^ - Ry/Vri^ , (B6) 

with the requirement that Rh = 1 (thus ji = k 2 ) in order to have the induced metric g^ v \R H = r/^. The coordinate 
transformation can be written locally as 

R= l + jk 2 r 2 + 0(r 4 ) . (B7) 
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Therefore, the above quantities are 

L= h 2 r + 0(r 3 ) , (B8) 

Kv» = \k 2 r^ v + 0(r 2 ) , (B9) 

from which we confirm that (3 2 = and verify also the correctness of (3.29). This is because i^„| r= o = 
(— N 2 , A 2 , A 2 ,A 2 ), with N 2 — A 2 because of the Minkowski symmetry. Then, the relation (3.29) gives 

A 2 = N 2 = . (BIO) 

Moreover, from (B8) we obtain (3, in the case that 9 is normalized so that 9 € [0,27r). 
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